Statistics of Transmission Eigenvalues for a Disordered Quantum Point Contact 



G. Campagnano, O. N. Jouravlev, Ya. M. Blanter, and Yu. 
Department of NanoScience and DIMES, Delft University of Technology, Lorentzweg 1, 

(Dated: February 2, 2008) 



V. Nazarov 

2628 CJ Delft, The Netherlands 



o 
o 



I 

O 

o 



> 

00 

o 
o 
o 



I 

O 

o 



We study the distribution of transmission eigenvalues of a quantum point contact with nearby 
impurities. In the semi-classical case (the chemical potential lies at the conductance plateau) we find 
that the transmission properties of this system are obtained from the ensemble of Gaussian random 
reflection matrices. The distribution only depends on the number of open transport channels and 
the average reflection eigenvalue and crosses over from the Poissonian for one open channel to the 
form predicted by the circuit theory in the limit of large number of open channels. 



PACS numbers: 73.63.Rt, 73.23. Ad 



I. INTRODUCTION 

A quantum point contact (QPC) is one of the refer- 
ence systems of mesoscopic physics. The experimental 
discovery of conductance quantizaton^ triggered further 
research which contributed much to our modern under- 
standing of nanoscience. QPC is a constriction defined in 
a 2DEG by gates. The width of the constriction can be 
changed by the voltage applied to these gates. In the adi- 
abatic regime, if the distance between the gates changes 
slowly compared to the wavelength of an electron, the 
theoretical description is readily obtained?^?-. The two- 
dimensional motion of an electron confined between the 
gates is equivalent to one- dimensional scattering of an 
electron at a potential barrier. The height of the bar- 
rier is different for different transport channels. Semi- 
classically, the electron is fully transmitted if its energy 
exceeds the top of the barrier in a given channel, and is 
fully reflected otherwise. Thus, semi-classical transmis- 
sion eigenvalues of a QPC are strongly degenerate: One 
has a finite number of transmission eigenvalues equal to 
one and an infinite number of transmission eigenvalues 
equal to zero. This picture would manifest experimen- 
tally in the precise quantization of conductance as a func- 
tion of gate voltage. 

In real experiments, this degeneracy is lifted. A brief 
glance at any of many available experimental studies 
shows that conductance does not rise in ideal steps. The 
question whether the transmission eigenvalues are degen- 
erate is also important for a number of other reasons. For 
instance, if a QPC is prepared in a superconducting ma- 
terial, discrete subgap (Andreev) states develop^. These 
states describe quasiparticles localized around QPC. The 
number of these states equals the number of transport 
channels, and their energies are expressed via transmis- 
sion eigenvalues T„, i?„ — 1 — Tn sin^ Lp/2, with A 
and being the superconducting gap and the phase dif- 
ference across the QPC. If the transmission eigenvalues 
are degenerate, the Andreev levels are also degenerate. 
Thus, any small perturbation would lift this degeneracy 
and produce a number of states with very close energies. 
Such a perturbation would then drastically affect prop- 
erties of the system. 



An obvious candidate for this degeneracy lifting is 
quantum tunneling across the top of the barrier. Indeed, 
for a given energy there is a range of gate voltages when 
one transport channel has a transmission eigenvalue be- 
tween zero and one - a partially open channel. All other 
transmission eigenvalues are also modified by the quan- 
tum tunneling: They get an exponentially small correc- 
tion. Thus, quantum tunneling leads to the rounding of 
the conductance steps as a function of gate voltage, but 
only provides exponentially small splitting of Andreev 
states. 

In this paper, we study how the degeneracy of trans- 
mission eigenvalues is lifted by the scattering on impuri- 
ties, which are always present in and around the QPC. 
Properties of a disordered QPC have been investigated 
(see Refs. HHSIHibj mostly in relation to the disorder 
smearing of conductance steps or evolution of conduc- 
tance fluctuations in ballistic regime. In contrast to the 
previous literature, we investigate the case when the con- 
ductance of the QPC is only slightly modified by the im- 
purities, or, in other words, the impurity-related splitting 
of transmission eigenvalues is much less than one. This 
regime is realized for low concentration of impurities. In 
this situation we can disregard quantum effects like reso- 
nant tunneling through impurity states or Kondo effect. 

Our main result is that in this regime, reflection am- 
plitudes are Gaussian distributed with zero average and 
second-order correlation function which does not depend 
on the channel index. This provides us with a new class 
of random matrix theory. The results for the distribution 
function of transmission eigenvalues are universal — they 
only depend on the number of transport channels and on 
the average reflection eigenvalue. All other information 
can be extracted from these two parameters. 

The paper is organized in the following way. In Section 
mwe treat a disordered QPC in the adiabatic approxi- 
mation. In Section llTTl we introduce the scattering matrix 
and show that in the expansion up to the second order 
in disorder potential closed channels do not contribute to 
the properties of transmission eigenvalues of open chan- 
nels. Section ITVI finalizes the quantum-mechanical calcu- 
lation of reflection coefficient and conductance of a dis- 
ordered QPC. We then turn to the classical (Boltzmann 
equation) consideration, which facilitates the considera- 
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tion of the diffusive regime (Section |VJ. 

In Section rVn we discuss noise properties of disordered 
QPC. Finally, Section rVIII is devoted to the distribution 
function of transmission eigenvalues. For one open trans- 
port channel, we calculate this distribution function an- 
alytically by performing the disorder averaging directly. 
In the limit of large number of open channels, we obtain 
the distribution function by means of the circuit theory^*', 
which presents a disordered QPC as a pure QPC and a 
diffusive resistor connected in series. For intermediate 
numbers of open channels, we perform a numerical sim- 
ulation based on random matrix theory. 



II. MODEL OF QPC WITH IMPURITIES 

We describe the QPC as a constriction between two in- 
finitely high walls^ separated by the distance d{x) (Figure 
1). A more physical model would take into account that 
the transverse profile is not sharp^. Since in this paper we 
employ semi-classical approximation (do not discuss the 
rounding of conduction steps), the results do not depend 
on the details of the potential profile. For this reason, we 
use the simpler model. The Schrodinger equation. 



^p{x,y) ^ Eip{x,y), (1) 



is supplemented by the boundary conditions, 
^{x,y^±dix)/2) = 0. 

Here 

V{x, y) = ^ v{x -Xi,y- yi), 

i 

with V being the single impurity potential, and the sum 
is taken over impurity positions. 

If the width of the constriction d{x) changes smoothly, 
we can employ the adiabatic approximation and separate 
the transverse motion. 



The transverse wave functions ip^i^ (y) that satisfy the 
boundary conditions are 



d{x) 



d{x) 



y 



d{x) 



Substituting this into Eq. ^ and disregarding the terms 
containing the derivatives of d{x), we obtain a one- 
dimensional equation for the longitudinal wave function. 



d^ 



E 



i)n{x) = -^Km(a;)0m(x), 

m 

(2) 



with the channel-dependent effective potential barrier 

'"^•''^ ~ 2md^{xy 
and the matrix element of the disorder potential. 



^nm ('^) 



d{x)/2 



-d(x)/2 



dy^'^:\y)V{x,y)v^-^\y). 



Eq. @ is the generalization of the equations previously 
written in Refs. ,2^ to the case of disordered QPC. 

In the semi-classical (WKB) approximation, in the ab- 
sence of disorder, for each transport channel, electrons 
with the energies above (below) the top of the barrier 
are perfectly transmitted (reflected). This approxima- 
tion breaks down if the energy of an electron coincides 
with the top of the barrier. In this paper, we do not 
consider this case. The wave function of an ideally trans- 
mitted electron is 



/Pn(oo) 



■ exp 



Pn{z)dz 



(3) 



Pn{x) 

with the channel-dependent momentum 
p„(x) = [2m{E - €^{x))Y'^ 



III. CORRELATORS OF THE SCATTERING 
MATRIX ELEMENTS 

We proceed by introducing the scattering matrix, 

S-( ^ ^ 

which is unitary, S — 1, due to the current conserva- 
tion requirement. At zero temperature, conductance of 
the system is expressed via Landauer formula. 



G = GgTrFt = a 



Tn 



where Tn (of interest in this paper) are the eigenvalues 
of the matrix rt, and Gq = /nh is the conductance 
quantum. Without impurities, the matrix tH is diagonal, 
with the elements describing the transmission of an elec- 
tron in the same open transport channel equal one and 
all others equal zero. In this case, the conductance is 
Go = GqN, with N being the number of open transport 
channels. 

To treat the effect of disorder, it is more convenient 
to investigate the matrix f^f, with the eigenvalues {re- 
flection eigenvalues) Rn = 1 — Tn- In the following, we 
calculate the correction to the transmission eigenvalue 
Rn due to disorder. For this purpose, we consider the 
perturbation expansion of the reflection matrix f up to 
the second order in the disorder potential. 



s(o) 



?(2) 
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Let us now separate open and closed channels, 



oo ' oc 



where the submatrices rap, a, (3 = o, c describe reflec- 
tion from/to channels of different type (o and c stand for 
open and closed). The reflection eigenvalues are found 
from the secular equation, det (f^f — — 0, or, equiv- 
alently. 



det( 

' nr' OO ~ ' rr' C 



Rl '^oo'^oc + '^to'^cc 
3 '^oc^oc "i" '^cc^cc ^ 



: 0. 

(4) 

We now expand this equation in powers of the disorder 
potential. For open channels, the reflection eigenvalues 
R are expected to be of the second order in disorder. We 
use now the identity 



det 



A B 
C D 



dei{A- BD-^C) detL>, 



and expand the first determinant taking into account that 
f'^oc = fco = Tec'' = 0. The terms of zeroth order in the 
disorder potential cancel, since f'i^'^^fi^^ = 1. Terms of 
the first order do not appear, and in the second order one 
has det{fi])^ ri]) — Rl) = 0. Thus, closed channels have 
no effect on transmission of open channels. In the rest 
of the paper, we drop the subscript oo and operate only 
with the quantities related to open channels. 

Elements of the matrix f are random quantities. In the 
next Section, we characterize their statistical properties. 
We show that they are Gaussian distributed with zero 
average. Thus, it is enough to specify the correlation 
functions 



' ij ' kh I 



where the average is performed with respect to disorder. 
We show that the only non-zero correlator is ^kij"*!^^! 
which is the probability for an electron coming in the 
open channel i to be scattered to the open channel j. 
As a matter of fact, the result does not depend on the 
indices i and j. We then relate this reflection probability 
with the correction to the conductance. 

Thus, the set of eigenvalues i?„ describing open chan- 
nels is obtained by diagonalizing the finite-size random 
matrix f^f. This matrix is Gaussian, and depends only 
on one parameter — average reflection coefficient. This 
novel random matrix theory in later Sections provides a 
novel distribution of transmission eigenvalues. 



IV. SCATTERING MATRIX APPROACH 

To calculate the correction to the reflection amplitudes 
we consider the Green's functions of Eq. Q , 



2m dx^ 



Gn{x, x') = -5{x - x'). 
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FIG. 1: (a) Layout of a disordered QPC. (b) An equivalent 
circuit representing a disordered QPC as a clean QPC and a 
diffusive resistor in series. 



Solving this equation, we find the Green's functions for 
open channels. 



Qn{x,x') = 



ih 



: exp 



dzpn{z) 



m^Jpn{x)pn{x') 

The formal solution of Eq. (0) takes the form 

) = + j dx'g 



(5) 



(6) 

with 4)n^ being the solutions in the absence of disorder 
(Eq. m with the zero right-hand side.) In the first order 
in y, we obtain 



(j^nix) ^ 4^\x) / dx'gn{x,x')V,,„,{x'^t\x'). 

Substituting Eqs. © and © into Eq. Q, we find 

oo m^pn{x')pm{x') 



(7) 



X exp 



dz {pn{z) +Pyn{z)) 



VnM)- (8) 



For the Gaussian distribution of disorder, the reflection 
amplitudes r„m are also Gaussian distributed. This dis- 
tribution is fully characterized by the pair correlation 
function. If the impurities are not located in the con- 
striction, the momenta Pn{x) and pm{x) in Eq. (jH)) can 
be replaced by their values taken at x ^ cxd, which is 
independently of the channel index. The impurity aver- 
aging is straightforward. As anticipated, all the averages 
of the type {f*jfkh) turn to zero due to the oscillating 
behavior, except for the term {\r„ 



I nm l) , 



where w(7r) is the Fourier transform of the single impurity 
potential with the momentum transfer 2kp^ and is the 
concentration of impurities per unit area. 
The correction to the conductance reads 



N 



{5G) = -Gq J2 (I^ 




and in the case of large number of open channels N can 
be written as 



(SO) 



(10) 



For further reference, we identify the average reflection 
eigenvalue (i?) by means of Landauer formula, {5G) = 
-GqN{R), 

h vf, d 

The correlation function of the reflection amplitudes can 
then be expressed via only one parameter (i?) , 



(|rn«H = ^(i?> (1 + %^ 



(11) 



V. BOLTZMANN EQUATION 

In this Section we analyse transport properties of a dis- 
ordered QPC in the framework of the Boltzmann equa- 
tion. This approach allows for an extension of the anal- 
ysis to the diffusive case, when the mean free path be- 
comes much smaller then the length of the system. In 
this case the second order perturbation expansion breaks 
down, and the treatment of previous Sections can not be 
applied any more. 

We model the system as a two-dimensional disordered 
wire between ideal reservoirs. To take into account the 
constriction, we recall that without disorder, in the quan- 
tum treatment only the channels with low index are ide- 
ally transmitting (open). In the language of classical 
physics, these channels correspond to electronic modes 
propagating with small incident angle a (Figure 2). To 
implement this feature in our model, we assume that all 
electrons propagating with angles smaller (larger) then 
ao ^ 1 are perfectly transmitted through the cross- 
section X = (to model the position of QPC). Electrons 
with higher angles are reflected from this cross-section. 
The Boltzmann equation reads 



t;.V,/(r,p) =/[/], 



(12) 



where f{r,p) is the distribution function of electrons, 
and /[/] is the collision integralii. 



m 



27m,- 



d^p' 



■[f{r,p')-f{r,p) 



FIG. 2: Disordered QPC modeled classically. Electrons with 
the angle a lower than the critical angle ao are transmitted 
through the cross-section a; = 0, others are reflected from this 
cross- section. 



x\i,ip~p')\^Sieip')-e{p)). 

Since only electrons with energies betwen Ep and Ep + 
eV, with V being the applied voltage, contribute to the 
net current, the absolute value of the momentum p is 
fixed to lie at the Fermi surface. We are only interested 
in the angular dependence of the distribution function, 
f{x,a). The Boltzmann equation (|12|) is supplemented 
by the boundary conditions. 



/(x, a) = 1, cos a > 0, a; < 
/(x, a) — 0, cos a < 0, x > L, 



which state that electrons coming from the reservoirs 
are in thermal equilibrium. To take into account reflec- 
tion of electrons at the QPC, we introduce the further 
boundary condition, 

/(O, a) = /(O, TT — a), I cosal < cosao. 



The distribution function does not depend on the 
transverse coordinate y, and we thus rewrite Eq. H12|l 



vp cos a dxf{x, a) 



(13) 



Him 
2^ 



da' cos a' \v{a — a')\ [f(x, a) — f{x, a')] . 



We solve now Eq. H13|l in the two limiting cases of low 
and high impurity concentration. The first case corre- 
sponds to the consideration of previous Sections based 
on the scattering approach. The second one, instead, 
allows us to investigate the diffusive regime. 



A. Non-diffusive regime 

For low impurity concentration, we replace the func- 
tion f{x,a') in the rhs of Eq. IjlSI) with the distribution 
function of the pure system. The remaining differential 
equation is readily solved, yielding for the distribution 
function at a; = and with the direction of propagation 



5 



cosa| > cosao, 



VI. NOISE 



/(0,a) 



1, cos a > 

aoLnim\v{a)\'^ 



vpirTi cos a 



cos a < 



Now we calculate the current density, j = 
{ev{Ep)vp /2ti) ^ dacosaf{Q,Q), with v{Ef) being 
the density of states at the Fermi energy. Only the 
angles |cosa| > cosao contribute, and we obtain the 
conductance. 



2ao — 2a, 



2 riiTn 



2Trh vp 



d, (14) 



where we have taken into account that ao <C 1. The 
second term represents the correction to the conductance 
due to impurity scattering. Relating the number of open 
channels to the critical angle ao, ao = 2NTr/kpd, we 
find that this correction is identical to Eq. (|10|l . which 
we obtained from the scattering matrix approach in the 
limiting case iV 3> 1. 



B. Diffusive regime 



In a two-terminal system, zero-temperature shot 
noisei^ can be expressed in terms of reflection eigenvalues 
in the following way. 



S 



i?ri(l — Rn) 



In a clean QPC, all transmission eigenvalues in the semi- 
classical regime are either zero or one, and the current 
is noiseless. Quantum tunneling away from the conduc- 
tance steps only brings exponentially small contribution. 
Thus, we expect a drastic effect of disorder on the shot 
noise of a QPC. Indeed, up to the second order in the 
disorder potential, one writes 



S 



2e^\V\ 



-2e\V\SG = 2eGQN\V\{R). 



(16) 

Eq. (|16|l corresponds to the notion of Poissonian stream 
of reflected particles, which is, indeed, expected in the 
case of good transmission. 



For high impurity concentration, we introduce the 
transport relaxation time r, defined asii 



-3 . (ia|5(a)| (1 — cos a). 
2Trh Jo 



In the linear regime, the collision integral takes the form 
/[/] = —{f — f)/T, where / is the distribution function 
averaged over the angles. 



1 



2ir 



27r Jo 

Solution of the Boltzmann equation in this case is given 
by 



rl \ f ''^^ \ f X I 



with I — vpT being the mean free path, I <C L. Calculat- 
ing the current, we find that the conductance is 



_ e^vpv{Ep)d I I L\ ^ 



G , . , 

27r \2ao ttI 

Identifying the conductance of the diffusive region, 

e'^vpv{Ep)ld 



(15) 



Gd — 



2L 



and of the clean QPC, 

e'^vpv{Ep)d 



G 



QPC 



2-K 



2ao = GqN, 



we see that Eq. (|15|l represents a series addition of the 
conductances of the clean QPC and a diffusive resistor. 



VII. DISTRIBUTION FUNCTION 

Now, we turn to the calculation of the distribution 
function of transmission eigenvalues, 

where the sum is taken over all open channels, and dis- 
order averaging is performed. The distribution function 
is normalized so that its integral is the number of open 
channels N. For a clean QPC, p(T) = N6{1 - T). 

We first consider the case of one open channel and 
perform the disorder averaging directly. Then, we obtain 
analytical results in the limiting case of very large number 
of channels based on the circuit theory: A disordered 
QPC is presented as a clean QPC connected in series 
with a diffusive resistor. For intermediate values of N , 
we were not able to obtain analytical results. Instead, 
we perform a simulation based on the notion of random 
reflection matrices. 



A. One open channel 



For one open channel, we perform a "brute force" dis- 
order averaging. The reflection probability is R — |rp, 
where we have suppressed the channel indices. In its 
turn, the reflection amplitude r is related to the poten- 
tial V by Eq. 0, with to = n = 1. It is a complex 
quantity, and we first calculate the distribution function 
of the transmission amplitude, defined as 

V{r) = {5 (Re(r - r[V])) 5 (Im(r - r[V]))) . 
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Representing delta-functions as integrals, we write ex- we obtain 
plicitly 



nr) 



duJiduj2 



, (2^)2 

X exp {iujiRe [r[V] - r] + itj2lni [r[V] - r]}) . (17) 

Performing the averaging with the Gaussian distribution 
and introducing the short-hand notation 

r[V] = / drF{r)V{r), 



V{r) 



dLJiduJ2 

^exp 



. (2^) 
iWiRer — ici;2liiir} 



i j dridr2 [tuiReFin) + cj2lmF{ri)] {V{ri)Vir2)) [u;iReFir2) + uj2lmF{r2)] 



r 



Now we use Eq. (jSJ and disregard the terms containing 
rapidly oscillating functions. Calculating the integrals 
over duji and c?W2, we finally obtain 



Vir) 



1 



exp 



(18) 



Since dR = 2\r\d\r\ = tt^M Rer d Imr, we can rewrite 
Eq. (|18l) in terms of the reflection eigenvalue, 



p(i?) = — exp 



R 



p(T) = p(l-i?). (19) 



Thus, for one channel the reflection eigenvalue is Pois- 
son distributed. This result is actually not surprising. 
Indeed, both real and imaginary parts of r are linearly 
related to the potential V. This means they are both 
Gaussian with zero average. Moreover, from Eq. Q it 
follows that they have the same dispersion, and we ar- 
rive to the Gaussian distribution for r (|18l) and Poisson 
distribution for R p9|l . These distributions are univer- 
sal: All the information about the type and amplitude 
of disorder is encoded in only one number, which is the 
average reflection eigenvalue (R). 



B. Circuit Theory 

If the number N of open channels is large (G ^ Gq), 
we can calculate the distribution function p{T) analyt- 
ically by means of the circuit theory developed by one 
of the authors-^^°. To this purpose, we represent a disor- 
dered QPC as a clean QPC connected in series to a dif- 
fusive conductor (Figure 1). Such a point of view, to our 
knowledge, was first adopted in Ref. for investigation 
of conductance fluctuations. The input parameters are 
the conductances of both circuit elements (connectors), 



Gqpc and Go- Each connector is subject to a phase dif- 
ference (j), which generates the pseudo-current /(0). The 
relation I{4>) is determined by the distribution function 
of transmission eigenvalues. 



1 



1 



/ I TT + 2iarccosh 



1 

Vt 



p(T) = ; 

In our case, the current-phase relations are 

• Diffusive conductor: /(0) = Gd'P 

• Quantum point contact I{(f>) = 2Gqpc tan(0/2) 

The circuit theory shows how these two elements can be 
combined to get the distribution function of the entire 
circuit. To this purpose, we introduce the phases: (f> in 
the left reservoir; zero in the right reservoir, and (to 
be calculated) in the point (node) separating the QPC 
and diffusive conductor. Thus, the QPC is subject to 
the phase difference (j) — 0, and the diffusive conductor 
to the phase difference 9. The pseudo-current must be 
conserved, from which we get the following equation for 



Gd0^2Gqpc tan((0-0)/2). 



(20) 



After solving this equation, we find the pseudo-current 
7(0) = G_D0(0), and eventually the distribution function. 

Eq. 1)20(1 can not be solved analytically for an arbitrary 
relation between Gd and Gqpc, and we restrict ourselves 
to the case of low impurity concentration, Gd ^ Gqpc- 
We obtain 



/(0) = Gd 



and the distribution function follows 
Gd 1 




P(T) 



27rGQ tVT^ 



(21) 
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for T > Tc, Tc = 1 - AGqpc/Gd, and zero for 
T < Tc- For further comparison with the simulation 
results, we calculate the average reflection coefficient 
(R) = Gqpc /Gd <C 1 and rewrite the distribution fmic- 
tion |(5lJ as 



p{R) 



N 



1 



27r(i?) 1 - i? V R 



- 1. 



(22) 



We see that in the limiting case N ^ 1 only channels 
with the transmission close to perfect exist: Reflection 
eigenvalues can only be lower than 4(i?), which is a small 
number. This is in contrast to the one-channel case, 
where all values of the reflection coefficient are permit- 
ted. To demonstrate the crossover between these two 
limiting cases, we perform a numerical simulation for the 
intermediate number of open transport channels. 



A 

on 

V 
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V 

2 




FIG. 3: The distribution function of reflection eigenvalues 
p{R) for 1, 5, 20 and 50 open channels. It is calculated using 
an ensamble of 10^ random matrices. The dashed line is the 
expression 1221 valid in the limit N oo. 



C. Intermediate case 

For any number of open channels, the matrix elements 
of the reflection matrix, r„m are random quantities with 



Gaussian distribution. They are zero on average, uncor- 
related, and characterized by the dispersion Eq. 
© shows that this dispersion is with the factor of 3/2 
greater for diagonal elements that for non-diagonal ele- 
ments, and otherwise does not depend on n and m. Then, 
the problem has only two parameters — the number of 
open transport channels N, and the dispersion of matrix 
elements of the reflection matrix. The latter parameter 
is related to the average reflection eigenvalue. 

Thus, we model the reflection matrix as a Gaussian 
random matrix with independent elements Both 
real and imaginary part of r„„i are taken to be Gaussian 
distributed with the same dispersion (independent on n 
and to). The simulation involves 10^ matrices. Results 
for 1, 5, 20, and 50 channels are plotted in Figure 3. 
For better comparison with Eq. (|22|l . we multiply the 
distribution function with i?^/^, so that it tends to a 
constant value at i? — *■ 0. For one channel, the simulation 
result perfectly reproduces the Poisson distribution l(T^ . 
whereas for 50 channels it is in a good agreement with Eq. 
(|22|l . For 5 and 20 channels, the simulations reveal an 
intermediate picture, with a tail at large R. Oscillations 
in the distribution function, pronounced for 5, 20, and 50 
channels, are related to the Wigner-Dyson correlations of 
eigenvalues of Gaussian random matrices. 



VIII. CONCLUSIONS 

The main result of our paper is the distribution func- 
tion of reflection eigenvalues of a disordered QPC. We 
have shown that for one open channel, one obtains Pois- 
son distribution, and in the limit of infinitely many chan- 
nels only very small reflection eigenvalues are allowed. 
We also performed numerical simulations for the inter- 
mediate regime. In all the cases, the results are universal 
in the sense that they only depend on the number of open 
transport channels and on the average reflection coeffi- 
cient (R). To calculate (i?), we developed a quantum- 
mechanical theory, and a classical calculation based on 
the Boltzmann equation. For low impurity concentra- 
tion, classical and quantum-mechanical results are iden- 
tical; additionally, classically one can treat the diffusive 
regime, to find that the conductance of a disordered QPC 
is given as a series resistance addition. We stress that al- 
though the results for (R) are model dependent, as care- 
fully investigated previously by Glazman and Jonson^, 
the form of the distribution function only uses (i?) as a 
parameter, but does not depend on specific model any 
more. It can be regarded as a result of the novel random 
matrix theory with Gaussian reflection matrices. 

This work was supported by the Netherlands Founda- 
tion for Fundamental Research on Matter (FOM). 
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